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Abstract. We have discovered three non-power infinite series representations for Bessel 
functions of the first kind of integer orders and real arguments. These series contain only 
elementary functions and are remarkably simple. Each series was derived as a Fourier 
series of a certain function that contains Bessel function. The series contain parameter b 
[ by setting which to specific values one can change specific form of series. Truncated series 

retain qualitatively behaviour of Bessel functions at large x: they have sine-like shape with 
decreasing amplitude. Derived series allow to obtain new series expansions for trigonometric 
functions. 



Introduction 

Bessel functions are widely used within mathematical physics, which explains everlast- 
ing interest to these functions both in physics and mathematics. Key properties of Bessel 
functions were articulated about a hundred years ago. Those results are mainly short, good 
looking and represent classic results of special function theory. They are collected in numer- 
ous treatises (e.g., see [T]). 

Often works on Bessel functions stem from the needs of physics and are focused on the 
properties of Bessel functions required to deal with a particular physical problem. That 
was the reason to study Bessel functions of large order and argument [2H1] which arise in 
astrophysics [51E]- 

Most works on Bessel functions — discussing algorithms for Bessel functions calculation [21 
IThTU]. identities containing Bessel functions [TlTfH] and others — yield cumbersome results. 
The result of present work, obtained for Bessel functions, is concise and simple. 
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1. Infinite series representations for Bessel functions of real order 
We start from three known identities [T5] 2.12.21]): 

j\" +1 J u (bx)sm{yVT^)dx= ^yb v {^¥^Y^ J v+ ,{^/¥^) , (la) 



(lb) 



r 1 J v {bx) r tt ( jWTtf-y \ f V^T7 + y \ 

f ^=^ C0S ( yVY ^~ 2 ) dx = ^(V^^T^-J^iVb^), (1c) 

where b > 0, y is any real number and v > — 1. Left sides of Eqs. ([[]) can be transformed 
into another form by introducing a new variable t = \/l — x 2 : 

/ t(VT^¥y 'j v {bVl^) sin (yt) dt = y/2^yb v (y/V> + y*) " 5 J u+ i (y/p + y 1 ), (2a) 
cos(yt) dt = 7i J u/2 J u/2 - , (2b) 



i 



J u (bVT^¥)(Vl :r Py cos(yt) dt = yfhtb v {>JW + y 2 ) " * J y +± {\/b 2 + y 2 ). (2c) 

Here we extended interval of integration over t from [0, 1] to [—1,1]. 
Eq. (j2j) can be presented in general form 



(b, f )e»* dt = / W (6, y) , a = A, B, C, (3) 



where 



F^\b,t) = -^(vr^yjuibvr^), 

f< A \b,y) = ^yW{^¥^ 2 )- u ~'h v M^¥^y~ 2 )., 



/J J (6, J/) = 7rJ^ /2 I J J u/2 I 

Ff\b,t) = j„(bVT^¥)(VT~p)\ 

fP(b,y) = v^(^^n~^ + i(^+7). 



(4a) 



(4b) 



(4c) 



If y — /ot, one can consider integrals on the left side of Eq. ([3]) to be calculations of the 
Fourier series coefficients of function Fj) a \b,t) defined on interval t G (—1, 1) and extended 
to a periodic function on BL Therefore one can present F^ (b, t) in form of the Fourier series 

FW(M) = - Y fi a) {bM)*-*«, I e hW l f ° = 1 ^ " 1 n < " - °' R n ( 5 ) 
2 ^ ^ v ' y t e [-1, 1] if a = A and z/ > 0, or a = B, C. w 

k=—oo 
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In case of a singularity at t = ±1 on the left side of Eq. (J5J) corresponding series on the right 
side will also give singularity. In case of the removable singularity at t = ±1 we have to 
remove singularity by redefining the function at these points and the equation will remain 
true. Series in Eq. fl5]) gives zero at points t = ±1 if a = A and — 1 < v < 0. All Fourier 
series (jSJ) are valid while F„ A ' B ' c \b,t) are differentiate with respect to t at t G (—1, 1). 



Coming back to variable x = yl — t 2 Eq. (JSJ) is transformed into 

1 oo 

Fl a) (b, VT^) = - f!) a) {b,k7c)e- lkwVI ^, where 



2 

k=— oo 



sc G (0, 1] if a = A and - 1 < v < 0, 

x G [0, 1] if a = A and z/ > 0, or a = B,C. 

Exchanging 6 and x we get the following series expansions at x > 



(6) 



1 oo 

(x, VT^) = - fi a \x,kn)e- lk7rVT ^ 1 , where 



2 

fc=— oo 



b G (0, 1] if a = A and - 1 < v < 0, 

b G [0, 1] if a = A and z/ > 0, or a = B, C, 

which in every particular case of a has the form 



(7) 



VVl^VUbx) = £ f?\x, kn) sin (WT^F), £ | fH 1 * " " °' (8) 

fe=i [ ' J ' 

which we will call (A)-case, 

^^ = E^/f ) (^^)cos(A;7rv / r^), &G[0,1], (9) 

k=0 

which we will call (B)-case, and 

oo 

b"J v {bx) = J2 £ kfP{x,k7r) cos (knVT^), &G[0,1], (10) 



k=0 



which we will call (C)-case, where 



Eqs. USD, ©, (HUD are valid at x > and v > -1. 

2. Infinite series representations for Bessel functions of integer order 

Derived series representations merit detailed consideration at integer v as in this case they 
get series representations through the elementary functions. Series representations for Bessel 
function of integer order are following: 

(A)-case 

oo 

b n J n (bx) = J2(^)gl A \b)ti A \x,kn), be (0,1), n>0, x > 0, (12) 
k=i 



(B)-casdi 

^ = E^f C) W/f ) (^^), be[0,l], n>l, *>0, (13) 



b 

(C)-case 



k=0 



b«J n (bx) = £ e^(6)/(°)(x, far), ^ M * > 0, (14) 

where 



fc=0 



^(&) = ^^[/^ and 3r\b) = cos (WT^). (15) 

We imply that (0, hir) = lim^o fn (x, kir) (a = A,B,C) if singularity at a; = is 
present. 

However, we need to perform some transformations in (B)-case and at n = 2m in order 
to get elementary functions on the right hand of Eq. (Q. Application of recurrence relations 
for the Bessel functions pQ 

4m j /, X J2m-l( bx ) , ^2m+l(te) . , ,., R s 

— J 2 m{ox) = 1 , m > 1 (16) 

trx 

to equation ([9l) results in representation^ 

^J2 m (bx)=^eJW(b)f$(x,kir), be (0,1), x>0, m > 1 (17a) 



fc=0 



/■(B) 
'2m 



(X,?/) = X • 7T 



a/x 2 + y 2 - y\ ( \Jx 2 + y 2 + y 

J(2m-l)/2 ~ J(2rn-l)/ 



Jx 2 +y 2 -y\ (Vx 2 + y 2 + y 

+ J(2m+l)/2 I ~ I J(2m+l)/2 



(17b) 



We can extend derived expansions to the region x < according to the parity of the Bessel 
functions of integer order: J n (— x) = (— l) n J n (x). 
With notation 

<p k = vV + (far) 2 (18) 
the series expansions for the Bessel functions of order 0, 1, 2 are presented below for every 
particular case. 



-^Here we will not receive series representation for Jo (a:) through elementary functions (see Eqs. © 
and (|4bl) ) . Therefore we consider only case of n > 1 . 



2 Left side of this equation should be taken as lim^o at b = 0. 

3 We should note that there are numerous ways to derive Bessel functions of even order through Bessel 
functions of odd order. For instance, recurrent relations [1] give another two equations: 

1 dJ 2m+ i(bx) 



J2m(bx) = (2m + 1) J 2m +i{bx) + 



b dx 



J2m{ox) = (2m - l)J 2m -i{ox) - - . 

o da; 
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(A)-case: 



J (bx) 
Ji(bx) 
J 2 {bx) 



k=l 



2 h 9k \b) 



k=l 

„2 00 



.•''"v^ , ^ (A),^Vk{<Pk 2 - 15) cos^ fc + 3(5 - 2y? fc 2 ) siny? fc 



fc=l 



(19a) 
(19b) 
(19c) 



where b G (0, 1). 
(B)-case: 



2h 

J x {bx) = -J2s k g { k B ' C \b)[(-l) k - cosif k ]. 



X 

b 2 



k=0 

oo 



Mbx) = -J2^9 { k B ' C \b) [(-l) k (2 + x 2 ) - 2( Vk sin Vk + cos^)] , 



k=0 



(20a) 
(20b) 



where b G [0,1]. 
(C)-case: 



J (bx)=2j2e k9 l B ' c \bf-^. 



k=0 

oo 



J l {bx) = 2t-Y j e k gf' C \b) 



{b,c) (K , sm (f k - (f k cos (p k 



k=0 
2 00 



fc=0 



(3 - ip k 2 ) sin - 3ip k cos v? fc 



(21a) 
(21b) 
(21c) 



where b G [0, 1] for J Q (bx), and 6 G (0, 1] for J n (bx), n > 1. 

One can derive following properties of the fc-th term in series ffl2l . ffT3|) . f|14n at k — > oo: 



(kn)gi A \b)fi A \x,kn) < (kn) £ A \x , kn) 



n, m = 0, oo, 



(A;7r)/ 2 ^(x,A ; 7r)~2(-l) fc+m+1 (-^) 



/ x \ 



2m 



(k^Am+iix, kir) ~ 2 



V/C7T/ 

/_i\fc+m+l 2m+l 

1 1— J [x 2 + (2m + 2) (2m + 3 

kix / 



(22a) 
(22b) 
(22c) 



gf' c \b)fi B \xM) < f { n B \xM) 



n = 1, oo, m = 0, oo, 



, k + m x ( x ^ 2m+2 2m+l ! r 2 , /o 

^ +m — — - — ^— [x 2 + 4m(2m + 1 

2 2m+i Vfcvr/ 2 2m+l L V 



(23a) 
(23b) 



2in 



/£U*,*7t)~(-1) 

for) ~(-l) 



fc+m+1 



fc+m+1 



/ x \ 2m + 2 (2m + l)(2m+l) r9 , Nl , . 

I — ) -4 ^--r x 2 + 2m 4m + 3 , 23c 

\knJ (4m + 3)(2(2m + 1))! L v ;J ' v ; 

(23d) 



2 2m VfcTr/ (4m + 3)(2(2m + 1)) 

x 2m+l f x 2m +2 (2m + 2)! 



2 2m Utt/ (2(2m + 2))!' 



x 



2m+2 



x v 2m+2 (2m + 2) 



2 2m Utt/ (2(2m + 2))l 



gl B ' C) WP(x,kn) < fP(x,kn) 



(kir) 2 \kirJ 



, n, m = 0, oo 
[x 2 + 2m(2m + 1 



/£ } +1 (x,A;7r)~2(-l) 



k+m+1 



1 / X \ 2m +! 



fc7T V/C7T 



(23e) 

(24a) 
(24b) 

(24c) 



Notation ~ b k means that lim^oo |^ = 1. Eqs. fl22|) . fl23|) . fl24|) allow to conclude absolute 
convergence and uniform convergence with respect to parameter b at b G (0, 1) for series 
f[T2"j) . f[T3"j) . and f[T3"j) except for series f[T2"j) at n = (that is fl 1 9aj) ) which is conditionally 
convergent and is not uniformly convergent with respect to b at b G (0, 1). 

Due to the uniform convergence of series (1T21) with respect to b G (0, 1) at n > 1 we can 
take a termwise limit b — >■ 1 and hence put 5^(6) = 1 at 6 — > 1 in (A)-case: 



^n(x) = ^(^/^(X, fcTT), n > 1, X > 0. 



(25) 



fe=i 



Series (I25p at n = diverges. 

The series representations for the Bessel functions of order 0, 1, 2 at b = 1 are presented 
below. (A)-case: 



r / x_ 2^/, y Vk&k 2 - 15) cosy^fc + 3(5 - 2ip k 2 ) sinyg fc 

"2^J — ^X ^[KTT) j . 



fc=l 



(B)-case: 



M x ) = - ^2 £ k[(-l) k ~ cosy? fc ], 

k=0 
j oo 

Ji{x) = — ^2 £ k[(-l) h (2 + x 2 ) - 2(ip k smip k + cos^fe)]. 



(26a) 
(26b) 

(27a) 
(27b) 



fc=0 



(C)-case: 



J Q (x) 
Ji(ar) 



2 5> 

k=0 

oc 

2x e fc 



oo 



k=0 

oo 



sm y? fc - (p k cos y9 fc 

<£>fc 



(3 - y? fc 2 ) sin if k - 3f k cos tp k 



k=0 



(28a) 
(28b) 
(28c) 



According to some simple considerations, we conclude that series f|T2l . ([TBI . ( 1T%|) at all 
permitted values of parameter 6, and series ff25l) converge uniformly with respect to 2 € [0, A] 
where A is arbitrary large real number. 

One can derive another series representations for Bessel functions by taking 6^1 and 
scaling variable x by 1/6. So (A)-case gives the following series representation for Jq(x) at 
b = V3/2: 



Jo(s) = 4 £(-!)" [(271-1)*] 



n=l 



V> n cos ^ - sin j/y 

^n 3 



where 



^!+[(2n-l)7r] : 



(29) 



3. Infinite series representations for sine and cosine functions 

Derived series expansions allow to present numerous similar expansions for sine and cosine 
functions. We will present some of them. 

Lets consider Eq. (I20al) which we rewrite in the following way: 



x- 



J^bx) 



2^£^f C) (6)[(-l) fe -cos^] : 



(30) 



k=0 



where <fk is presented in Eq. ( TTBl . Series ( 130|) is uniformly convergent with respect to 
b G (0, 1) therefore taking a limit b — > in Eq. (I3"U1) gives us the following result: 



cosx 



-1) ~cos ip k 



k=l 



Taking b = in Eq. (12 lap we derive series for sine function: 



1 - 



sm x 



x 



2 £ 



k=l 



-1) 



fc Sm (p k 



(31) 



(32) 



It can be also derived by differentiating Eq. ( 13T|) with respect to x. 

Subtracting Eq. (12 7aj) and Eq. ( 128bl) we get another series for sine function: 



sin a; 

x 



- 2 E 

fc=i 



^ (*7T) 2 



x smyjfc 



(33) 



which converges faster than the previous one. 



Equations ( l31~j) . (|32|) . (|33|) can be differentiated any times that will give plenty of other 
series representations for sine and cosine functions. 

Conclusion 

To sum up we outline key properties of the series representations f|T^]) . f lTSj) . (TH1) . and fTSol) 
for Bessel functions of the first kind of integer order: 

• series (fl~2]) at n > 1, ( !T3|) . (fl^|) are uniformly convergent with respect to parameter 
b G (0, 1) at all x; 

• series ( !T2|) . (fl3]) . ( !T4|) at all permitted values of b, and series fl25l) are uniformly 
convergent with respect to variable x G [0, A], where A is arbitrary large real value; 

• series ( TT2|) at n > 1, ( TT5|) . (JH]) at all permitted values of b, and series (T2oT) are 
absolutely convergent; series (filZj) at n = (that is (I19al) ) is conditionally convergent. 

As series ffT3j) . ffT4|) at 6 = 1 and (125]) become alternating and /c-th terms approach zero 
monotonically at sufficient large k and x G [0, A] (see Eqs. (122]) . (123]) . (124]) ). the bonds for the 
error terms of infinite sums approximations by truncated series to K-term are absolute value 
of (f^+l)-th term. As one can see, (i^+l)-th terms are monotonically decreasing functions of 
x in all three cases. Besides, truncated series (fl~2]) . f fT3|) . ffT4]) . ( 125]) have oscillating behaviour 
(sine-like shape) with decreasing amplitude at large x that is they behave exactly like Bessel 
functions. However, truncated series fail to provide correct asymptotic behaviour of Bessel 
function at x — > oo [T]: 
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